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FREE-MOLECULE  FLOW  IN  THE  AXIAL-FLOW  TURBO-VACUUM  PUMP 

by 

John  S.  Maulbetsch  and  Ascher  H.  Shapiro 

ABSTRACT 

The  rarefied  gas  flow  through  a  moving  cascade  of  flat  plates 
was  investigated  analytically  and  experimentally,  with  the  application 
in  mind  of  high- vacuum  pumping. 

The  theory  is  based  on  the  assumption  of  free-molecule  flow,  i.e., 
no  collisions  between  molecules.  The  essential  quantities  determining 
the  performance  of  a  blade  row  are  the  net  transmission  probabilities 
from  one  side  to  the  other.  These  transmission  probabilities  depend 
mainly  on  the  blade  angle  and  spacing-chord  ratio  and  on  the  ratio  of 
mechanical  speed  to  mean  molecular  speed.  In  this  report  the  transmission 
coefficients  are  calculated  by  numerical  solution  of  the  governing 
integral  equation®,  and  are  compared  with  the  results  of  earlier  Monte 
Carlo  calculations. 

Experiments  with  a  rotating  test  machine  having  either  a  single 
rotor  or  a  rotor  followed  by  a  stator  gave  results  in  good  agreement  with 
the  theory. 

A  special  means  was  developed  for  testing  stationary  cascades 
under  conditions  simulating  the  operation  of  a  moving  rotor.  This 
apparatus  was  tested  with  several  blade  geometries.  The  results  con¬ 
firmed  both  the  theory  and  the  feasibility  of  the  new  rarefied  gas 
wind-tunnel  concept. 
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NOMENCLATORE 

Aspect  ratio  (blade  span  to  width  of  passage  between  blades ) . 

Blade  chord. 

Average  molecular  speed  for  Maxwellian  distribution,  \)8KT/Tf 
Molecular  number  density. 

Effusion  flux,  number  of  incident  molecules  per  unit  area  per  unit  time. 
Pressure . 

Universal  gas  constant  divided  by  the  molecular  weight  of  the  gas. 

Blade  spacing  in  the  tangential  direction. 

Ratio  of  the  mean  mass  speed  of  the  gas,  relative  to  the  blades, 
to  the  most  probable  molecular  speed.  S  =  V,V 2RT . 

Temperature . 

w  Components  of  the  thermal  velocities  of  molecules. 

Magnitude  of  the  mean  mass  velocity  of  the  gas  relative  to  the 
blades. 

Blade  speed,  in  tangential  direction. 

Axial  component  of  mean  mass  throughflow  velocity  of  the  gas. 

Ho  Coefficient  (the  net  flux  of  molecules  through  the  cascade, 
expressed  as  a  fraction  of  the  flux  incident  on  the  upstream  side). 

Blade  angle  with  the  tangential  direction. 

Angle  which  the  mean  mass  throughflow  velocity  of  the  gas  makes  with 
the  axial  direction. 

Transmission  coefficient  (the  fraction  of  molecules  incident  on 
side  1  of  the  blade  row  which  ultimately  emerge  to  side  2). 


The  subscripts  1  and  2  refer  to  the  upstream  and  downstream  sides  of 
the  blade  row,  respectively. 
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Background  of  Research 

Until  recently,  it  had  always  been  assumed  that  methods  used  for 
pumping  or  compressing  gases  at  the  very  low  pressures  dealt  with  in 
vacuum  technology  would  necessarily  be  quite  different  from  those  normally 
employed  at  aerodynamic  pressures.  The  expected  pressure  ratios  per 
stage  associated  with  axial-flow  or  centrifugal  compressors  operating 
at  aerodynamic  pressures  are  so  lew  (1.05  to  1.2  for  the  axial-flow; 

1.5  to  2.0  for  the  centrifugal)  that  far  too  many  stages  would  apparently 
be  required  to  obtain  the  overall  pressure  ratios  required  in  high- 
vacuum  applications.  However,  at  very  low  density  levels  the  molecular 
structure  of  the  gas  may  no  longer  be  Ignored  in  considering  the  flow 
mechanism.  The  entire  operation  of  fluid  machinery  must  then  be  looked 
at  from  a  molecular  view  point.  The  normal  limitations  on  pressure 
ratio  —  adverse  pressure  gradients,  separation,  shock  waves,  surging, 
and  stalling  —  are  no  longer  meaningful.  That  large  pressure  ratios 
per  stage  can  be  achieved  by  axial-flow  machines  in  the  free  -molecule 
range  of  pressures  was  demonstrated  by  the  experiments  of  Becker  ^ 
and  of  Bablanian^2). 

The  present  project  was  initiated  under  sponsorship  of  the  Office 
of  Naval  Research  after  approximate  calculations  for  the  case  of  very 
high  blade  speeds  Indicated  the  possibility  of  extremely  high  pressure 
ratios  per  stage  and  gave  same  idea  of  the  Influence  of  the  blade  row 
geometry.  Our  object  has  been  to  provide  a  theoretical  and  experimental 
foundation  on  which  axial-flow  b laded  vacuum  pungis  may  be  designed  and 
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the  ir  performance  predicted.  Initial  results  vere  presented  in  our  first 

(3) 

report  to  the  Office  of  Naval  Research  under  this  contract.  These, 
together  with  some  additional  material,  have  been  summarized  in  two 
published  articles  ^ 

We  present  herewith  our  final  report  on  the  subject,  constituting 
a  definitive  account  of  the  performance  of  moving  cascades  of  flat  plates 
in  the  free-molecule  regime. 

Summary  of  Present  Investigation 

The  work  described  in  this  report  may  be  divided  into  three  sections. 

First,  the  theoretical  problem  has  been  solved  by  numerical  solution 
of  the  governing  integral  equation.  For  single  blade  revs,  this  made 
it  possible  to  extend  the  theory  over  a  wider  range  of  operating  parameters, 
and  also  to  eliminate  the  not  inconsiderable  statistical  errors  inherent 
in  the  Monte  Carlo  method  used  previously  in  reference  3. 

The  second  aspect  of  the  program  was  the  extension  and  refinement 
of  the  previous  experimental  Btudies,  using  the  existing  apparatus. 

The  third  and  most  novel  phase  of  the  research  was  the  testing  of 
free-molecule  cascades  in  an  adaptation  of  the  device  which  in  effect 
made  it  into  a  new  type  of  low-density  wind  tunnel.  This  allowed  us 
to  test  a  variety  of  blading  configurations  more  conveniently  and 
rapidly  than  would  be  possible  in  a  rotating  machine.  The  scheme  used  has 
application  to  a  wide  variety  of  lew-density  wind-tunnel  experiments. 
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ANALYTICAL  PROGRAM 


Analysis  for  a  Single-Blade  Rw 

The  Monte  Carlo  method  for  calculating  the  molecular  flews  through 

(3) 

the  machine  has  been  described  fully  by  Kruger  .  This  method,  essentially 
statistical,  was  chosen  because  at  that  time  it  was  thought  to  require 
less  machine  time  than  the  solution  of  the  integral  equations  in  terms 
of  which  the  problem  may  also  be  formulated.  This  was  undoubtedly  the 
case  for  the  multi-row  calculations,  for  which  the  molecular  velocity 
distributions  are  non-Maxwellian.  Furthermore,  the  direct  analogy 
between  the  individual  samples  of  the  Monte  Carlo  method  and  the  individual 
molecules  in  the  flow  made  the  method  look  particularly  attractive.  In 
fact,  for  problems  involving  a  complicated  geometry,  or  for  the  precise 
calculation  of  stages  in  series,  the  statistical  approach  may  well  be 
the  only  feasible  one. 

However,  in  some  cases  it  is  preferable  to  solve  the  integral 
equations.  For  single  rows  of  flat  plates,  such  as  are  considered  here, 
the  computer  time  involved  in  the  integral  solution  is  much  less  than  that 
needed  for  the  Monte  Carlo  solution.  Therefore,  in  order  to  cover  a  wide 
range  of  variables,  and  also  to  get  more  accurate  results  with  the  amount 
of  machine  time  available,  the  integral  equations  were  programmed  on  the 
IBM  7090  computer.  This  had  the  added  advantage  of  providing  certain 
intermediate  results  which  might  be  helpful  in  designing  better  blade 
shapes  for  such  varying  requirements  as  high  pressure  ratio  or  high 
incident-molecule  capture. 

Fundamental  Mode  of  Operation 


The  operating  characteristics  of  the  machine  may  be  understood  by 
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considering  two  regions  of  gas,  separated  by  a  row  of  flat-plate  blades 
(Fig.  1)  moving  with  the  speed  Vg.  We  shall  make  the  following  assumptions 

1.  The  upstream  and  downstream  regions  are  each  Maxwellian. 

2.  The  flow  is  free -molecular  with  respect  to  the  blade  dimensions. 
This  requires  that  the  mean-free  path  of  the  molecules  be  much 
greater  tnan  the  blade  dimensions,  and  allows  one  to  ignore 
Inter-molecular  collisions  within  the  blade  row. 

3.  The  temperatures  in  the  two  regions  are  equal.  This  is  Justified 
by  the  consideration  that  the  rate  at  which  energy  is  carried  to 
the  blades  and  the  casing  by  the  molecules  is  small  compared  to 
the  rate  at  which  it  can  be  conducted  or  radiated  away. 

4.  A  stationary  state  prevails. 

It  can  be  shown  from  molecular  flow  theory  ^  that  the  number  of 
molecules  incident  upon  the  blade  row  per  unit  area  per  unit  time  from 
a  Maxwellian  region  of  number  density  Tl  is  given  by 


N 


'n.'/IsT  I  l 


€  +  X 


(' +  Hi  *) 


(Eq.  la) 


where 


(Eq.  lb) 


and 


^  X  s 
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Now  let  2! | g  represent  the  probability  that  a  molecule  Incident 
upon  the  blade  row  from  region  1  will  ultimately  reach  region  2 
(perhaps  after  many  reflections  from  the  blades),  and  <Tg|  the  same 
probability  of  transmission  from  region  2  to  region  1.  Then,  for  zero 
net  throughflow,  the  number  of  molecules  passing  through  the  row  In 
opposite  directions  must  be  equal: 

H*,*  =  <*>• 2 

When  a  net  throughf low  exists ,  we  may  characterize  the  volume 
flow  by  means  of  the  Ho  Coefficient,  W,  defined  as  the  ratio  of  the 
actual  molecular  throughflow  to  the  wmMuimh  possible  throughf  low .  * 

This  maximum  Is  Nx,  and  therefore  the  requirement  of  conservation  of 
molecules  within  the  blade  rev  leads  to 


M,Zie  - 


(Eq.  3) 


The  effusion  flux  N  depends  on  according  to  Eq.  1,  and 

hence  depends  on  the  blade  speed  Vjj,  on  the  mean  mass  throughflow  velocity 
Vrji,  and  on  the  most  probable  thermal  speed,  V 2RT  •  However,  as  we  shall 
now  show,  t£c0S  is  in  fact  extremely  small.  From  the  geometry  of 

Fig.  1, 


(i  *  vT/v 


(Eq.  k) 


•The  Ho  Coefficient  Is  related  to  the  net  voluse  throughflow  by  the 
equation 


Net  Volume  Throughflow  «\j*T/*ir  ■  W 
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and,  by  definition, 


S«  V/[2RT 


(Eq.  5) 


Furthermore,  since  the  mean  mass  throughflcw  velocity  is  equal  to  the 
volume  flew  per  unit  area,  we  may  write 


>4-  ^ 


(Eq.  6) 


where  is  the  ratio  of  the  throughflcw  area  at  the  cascade  to  the 
throughflcw  area  a  mean- free -path  length  away  fran  the  cascade,  in  the 
region  from  which  the  molecules  are  incident  on  the  cascade.  Combining 
Eqs .4,5  and  6 ,  we  form 


m 

X  *  $  cjtS  A  as  — £L=.  W  K* 
'  VIV2RT  * 


(Eq.  7) 


Thus,  Eq.  la  may  now  be  written  as 


KaW  I 

‘X  *  — - —  -1-  g  +  oc 
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C  i  +  H  *)] 


(Eq.  8) 


Expanding  the  functions  on  the  right-hand  side  in  ascending  powers 
of  X  ,  and  simplifying,  we  get 


- 

■  tf ?l 


1+  l|?  X  +  X1 


4  T 

X  X 

T  +  3o 


(Eq.  9) 
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This  may  be  Inverted  to  give 


X  *  S  cet 


(Eq.  10) 


Nov  W  can  never  exceed  unity,  and,  in  any  practical  situation,  it 
is  likely  to  be  not  more  than  l/2.  In  our  experiments  with  the  one-row 
compressor  (single  rotor)  and  with  the  two- row  compressor  (single  rotor 
followed  by  single  stator),  was  equal  to  0.274.  Hence  3  cars  ^  was 
at  most  equal  to  about  0.04,  and  was  generally  far  less.  In  our  experiments 
with  the  stationary  cascade  apparatus ,  S  c^s  ^  was  always  less  than 
•035. 

The  foregoing  analysis  shows  that  for  most  purposes  it  is  sufficient 
to  assume  that  3  ces,  &  O  .  With  this  assumption  it  follows  from 
Eq.  1  that 


*,  w\j2'RT 

M  «  - ; — —  =  - 

2 '(IF  4 


(Sq.  11) 


Then,  since  the  temperatures  are  the  same  in  regions  1  and  2,  Eq.  3 
may  be  put  in  the  following  form  which  expresses  the  operating  per- 
formanance  of  the  row  insofar  as  it  relates  the  pressure  ratio 
across  the  row  to  the  volume  throughflov  parameter  W: 


nz 

W 

^21 

(Eq.  12) 
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Since,  as  shown  above,  the  net  throughflow  Imparts  only  a  very 
small  mean  mass  velocity  to  the  gas  In  the  axial  direction,  the  £  a. 
are  only  weak  functions  of  W.  Neglecting  this  dependence,  the  analytical 
problem  reduces  to  that  of  determining  the  values  of  2^  and  2j>(  as 
functions  of  the  dimensionless  blade  speed  3=  ^  /{Z%T  and  the 
blading  parameters  o L  and  s/b  .  Note  that  in  the  free-molecule  range, 
the  absolute  density  level  and  the  Knudsen  number  are  not  significant 
variables . 

Analytical  Formulation 

The  problem  may  be  divided  into  two  parts:  (1)  the  determination 
of  the  initial  incidence  distribution  of  molecules  from  the  Maxwellian 
regions  as  they  arrive  at  the  blade  surfaces;  and  (2)  the  evaluation  of 
the  ultimate  transmission  probabilities  for  molecules  within  the  blade 
row  after  they  have  been  emitted  from  one  of  the  blade  surfaces.  The 
analysis  is  illustrated  by  Fig.  2. 

The  fraction  0^  of  molecules  incident  upon  the  blade  row  from 
region  1  which  pass  directly  to  region  2  without  colliding  with  a 
blade  surface  can  be  computed  explicitly  in  terms  of  the  variables 

oL}  s/b>  ft  3. 

Next,  the  fraction  of  incident  molecules  which  impinge  directly 
upon  the  blade  area  dx  at  position  x  (or  dy  at  j  )  along  the 
surfaces  of  the  blades  may  be  calculated,  in  terms  of  the  so-called 
initial  Incidence  distribution.  These  are  characterized  by  the  density 
distribution  functions  7^  and  ,  such  that  djx  is  the  fraction 
of  all  molecules  entering  the  passage  from  region  1  which  impinge 
directly  on  the  area  dx  ,  and  similarly  for  "Vy  dtj 
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Now  let  to®  the  ultimate  probability  of  a  molecule's  reaching 
region  2  after  having  been  emitted  from  area  at  X  (after,  perhaps, 
a  number  of  collisions  vith  the  blade  surfaces),  and  let  be 

similarly  defined.  Then  the  overall  trazismisslon  coefficient,  g  5 
may  be  expressed  as 


(Eq.  13) 


where  it  is  to  be  understood  that  the  integrations  sure  over  the  entire 
length  of  each  blade. 

In  this  expression  all  quantities  may  be  calculated  directly  from 
the  known  quantities  except  for  the  functions  and  (7^  ($)• 

For  this  calculation,  consider  molecules  being  emitted  from  a  <kx  at 
X  on  the  left-hand  blade  surface.  Let  us  define  as  the 

fraction  which  passes  directly  from  Ax  to  region  2  without  collision, 
and  as  the  fraction  which  passes  directly  from  cLx  to 

some  dUj  at  ij  .  Then,  since  the  rate  of  molecular  incidence  on  any 
area  equals  the  rate  of  molecular  emission,  we  may  write 


^  =  + 


(Eq.  lUa) 


and,  similarly, 


V  *  ^2  +J 


(Eq.  14b) 
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In  this  pair  of  simultaneous  integral  equations,  the  functions 

^  and  F4*J  may  be  computed  from 
the  geometry  of  the  blading  and  the  cosine  law  for  diffuse  reflection. 
They  are  independent  of  the  speed  of  the  blade  row.  When  the  cosine 
law  is  valid,  it  may  be  further  shown  that  F x)  *  F  C*)  • 


For  the  special  case  of  flat  plate  blades,  with  which  we  are  dealing, 
certain  conditions  of  synmetry  prevail,  greatly  simplifying  the  problem. 
These  are  that,  at  equal  values  of  X  and  t j  , 


yi  ) 


(Eq.  15) 


Moreover,  these  four  values  of  0~  must  satisfy  the  conservation  relations 


^ +  F  ’ 


(Eq.  16) 


Using  Eqs.  15  and  16,  Eqs.  lUa  and  l4b  may  be  combined  and  reduced  to 

(Eq.  17) 


<y  -  (i-O-fv 


To  solve  this  equation,  the  integral  on  the  right-hand  side  was 
first  expressed  as  a  sum  of  twenty  terms,  based  on  the  trape¬ 
zoidal  rule  for  numerical  quadrature.  Then  the  function 
was  solved  for  by  iteration.  This  in  turn  gave  ^and  then 

Eq.  13  could  be  solved  for 

The  analytical  expressions  for  all  the  functions  defined  above 
are  to  be  found  in  Appendix  A.  The  functions  were  computed,  and  Eqs. 
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17  and  13  solved,  on  the  IBM  7090  computer  at  the  M.I.T.  Confutation 
Center,  using  numerical  quadrature  with  twenty  Intervals  for  each  integral. 


Results  of  Analytical  Program 

The  results  of  the  calculation  program  for  flat-plate  arrays  have 
been  collected  and  exhibited  in  several  ways: 

1.  In  Table  1  are  displayed  the  values  of  T  for  systematic  var¬ 
iations  of  the  geometric  parameters  oi.  and  $A>  and  of  the  blade- 
speed  parameter  S  .  Note  that  these  are  all  based  on  the 
assumption  S  ors  ^  =  0  .  Note  further  that  positive  values  of  £ 
correspond  to  Z ,,,,  while  negative  values  of  ^  correspond  to  . 

2.  Figures  3,  U,  5,  6,  and  7  show  the  same  information  graphically. 

3.  Figures  8,  9,  10,  11,  and  12  show  the  two  combinations  of  the  2  s 

which  are  indicative  of  the  performance  of  the  machine.  Note, 
from  Eq.  12,  that  8ives  the  pressure  ratio  for  zero 

throughflcw,  while  Z,^  gives  the  value  of  W  for  no 

pressure  rise. 

1*.  Figures  13,  lk,  and  15  are  cross-plots  of  Figures  8  through 
12,  showing  the  effect  of  spacing-chord  ratio  on  /^2| 
and  X|2  “  ^-2 1 

5.  Figures  16,  17,  and  18  are  similar  cross-plots,  showing  the 
effects  of  blade  angle  oL  . 

6.  Figures  19  and  20  show  certain  intermediate  information  regarding 
the  incident  flux  distributions,  lJy  and  V  and  the  ultimate 
transmission  probability,  <7^ ,  of  molecules  reflected  from  the 


blade  surfaces. 


In  examining  the  results  of  the  analysis,  one  must  do  so  in 

the  light  of  the  characteristics  of  operation  that  would  he  required 

of  the  individual  blade  revs  in  an  actual  compressor.  Equation  12 

shews  that,  for  high  pressure  ratio  per  stage,  the  ratio  *12  /  *21 

should  be  large,  while  for  large  pumping  capacity,  the  difference 

■jT|2.  -  £2|  should  be  large.  These  goals  are  generally  not  compatible. 

For  Instance,  the  highest  pressure  ratios  are  attained  with  *10° 

or  20  and  with  s/b«iA  rr  l/f.  However,  these  configurations  have 

low  values  of  the  throughflow  parameter  Zgi  • 

However,  practical  consideration  of  a  working  design  for  such  a 

compressor  indicates  that  the  situation  is  not  necessarily  an  unfortunate 

one.  In  considering  the  requirements  for  high-vacuum  technology,  one 

-8  -3 

normally  thinks  in  terms  of  absolute  pressure  levels  of  10  to  10  mm  Hg. 

or  lower.  If  a  compressor,  such  as  the  one  under  discussion,  were  to  be 

used  to  attain  such  a  pressure  level,  it  would  probably  be  operated  in 

conjunction  with  a  mechanical  forepump  with  a  blank-off  pressure  of, 

*•2 

at  best,  \0jA.  (10  ran  Hg.).  This  would  require  a  pressure  ratio  across 
the  machine  of  approximately  10^  or  more.  Such  a  pressure  ratio 
clearly  cannot  be  obtained  with  a  single  blade  row.  Operating  in  air, 

(with  a  most  probable  molecular  speed  at  68#F  of  approximately  1350  ft/sec.), 
the  best  blade-speed  ratio  one  could  hope  to  attain  is  about  unity.  Even 
at  zero  flow  this  would  require  at  least  six  or  seven  stages.  At  a 
zero-flow  condition,  the  most  efficient  design  would  clearly  be  one 
which  incorporated  identical  stages  all  providing  the  highest  possible 
pressure  ratio.  However,  since  even  at  "blank-off”  there  is  always  a 
net  throughflow,  and  since  one  must  pass  through  a  "pimp-down"  period, 
high  volume  flow  is  actually  as  important  for  the  machine  as  a  whole  as 
high  pressure  ratio. 
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Nov  the  mass  throughflcv  is  the  sane  for  all  stages  of  a  multi-stage 
machine,  and  thus  the  volume  throughflov  la  inversely  proportional  to  the 
pressure  level  upstream  of  each  row.  Hence,  for  a  multi-stage  machine, 
all  the  stages  need  not  perform  the  same  function.  The  stages  near  the 
outlet,  where  the  pressure  level  is  relatively  high,  should  provide 
high  pressure  ratios  per  stage .  for  in  this  region  the  Ho  coefficient 
is  nearly  zero.  Hear  the  low-pressure  end,  on  the  other  hand,  where 
the  volume  throughflcw  is  likely  to  be  quite  large,  blade  configurations 
should  be  chosen  to  accomodate  high  flow  and  still  provide  some  reasonable 
pressure  ratio.  Thus,  closely-spaced  blades  with  small  blade  angles 
are  desirable  near  the  high-pressure  end  of  the  machine,  while  more 
widely  spaced  blades  with  larger  angles  are  desirable  for  the  low- 
pressure  end  of  the  machine. 

Comparison  with  Results  of  Monte  Carlo  Calculations 

It  is  of  Interest  to  compare  the  conputed  results  presented  here 
with  the  corresponding  statistical  results  obtained  by  the  Monte  Carlo 
method  .  There  are  8l  cases  exactly  comparable  with  respect  to  the 
values  of  s/bj  c k*JL  3. 

The  two  sets  of  results  for  these  8l  cases  can  only  be  compared 
on  a  statistical  basis.  Accordingly,  the  percent  error  for  the  Monte 
Carlo  calculation  was  calculated  for  each  of  the  8l  cases,  on  the 
assumption  that  the  present  results  are,  in  fact,  correct.  These 
cases  were  then  arranged  in  ascending  order  of  error,  from  the  largest 
negative  error  to  the  largest  positive  error.  The  errors  were  then 
plotted  (Figure  21a)  in  the  form  of  the  integral  of  an  error  distribution 
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function:  i.e. ,  the  fraction  of  the  total  cases  with  an  error  less 
than  a  given  error  was  plotted  against  the  given  error.  The  smoothed 
integral  curve  of  Fig.  21a  was  then  differentiated  to  give  the  error 
distribution  function  of  Fig.  21b. 

On  the  assumption  that  the  result  of  a  Monte  Carlo  calculation  would 
be  normally  distributed  around  the  correct  answer,  the  number  of  samples 
used  vas  approximately  chosen  such  that  there  would  be  only  a  10.96$ 
probability  that  any  particular  result  would  be  in  error  by  more  than 
10$.  This  is  sufficient  information  to  construct  the  theoretical  error 
distribution  function  from  standard  probability  tables.  The  normal  dis¬ 
tribution  function  for  the  "expected"  standard  deviation  of  O'-  0.0625  is 
shown  in  Fig.  21b. 

The  actual  error  distribution  and  the  "expected"  distribution,  as 
exhibited  in  Fig.  21b,  are  in  general  agreement,  to  the  extent  of  giving 
credence  to  both  the  Monte  Carlo  results  and  to  the  present  results.  That 
value  of  the  fractional  error,  with  respect  to  which  half  the  errors  are 
greater  and  half  less,  is  about  -1/2$,  which  seems  quite  reasonable  for 
a  sample  of  8l  cases.  This  comparison  of  the  means  may  be  associated 
with  the  small  truncation  error  of  the  integral  calculations,  as  described 
belcw,  but  this  is  somewhat  speculative. 

The  important  conclusion  to  be  drawn  from  Figure  21  is  that  both  the 
Monte  Carlo  results  and  the  present  results,  reached  by  completely  different 
calculation  routes,  are  in  substantial  agreement.  The  implication  is  that 
both  sure  substantially  correct.  We  believe,  however,  that  the  present 
results  are  the  more  reliable  in  individual  cases. 

Truncation  Errors  in  numerical  Solution 

In  order  to  Investigate  whether  a  division  of  each  blade  surface 
into  twenty  tones  for  purposes  of  numerical  quadrature  of  the  Integrals 
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in  Eqs.  13  17  was  sufficient;  trial  calculations  were  also  for 

ten,  forty,  and  eighty  tones.  Frctn  these  trials  it  appears  that  the 
truncation  error  is  less  than  1/2$  . 


EXTENSION  OF  COMPRESSOR  RESEARCH 


The  second  phase  of  the  investigation  involved  the  refinement  of 
the  experimental  program  on  the  single-row  compressor  with  oc=  20*, 
s/b  =  0.46,  and  an  extension  to  a  two-row  machine  comprising  a  rotor  and 
stator  having  mirror- image  geometries. 

Single-Row  Compressor 

The  experimental  work  vas  performed  on  the  existing  apparatus. 

A  cross-sectional  layout  of  the  machine  is  shown  in  Fig.  22.  Figure  23 
is  a  schematic  of  the  associated  equipment  and  instrumentation. 

Although  a  description  of  the  apparatus  and  the  experimental  technique 

(3) 

appears  in  our  first  report  ,  it  is  repeated  here  for  the  convenience 
of  the  reader. 

The  test  section  consisted  of  a  cylindrical  steel  housing  in  which 

the  aluminum  compressor  rotor  was  driven  by  a  high-frequency  induction 

motor.  To  maintain  densities  sufficiently  low  for  free-molecule  flow 

in  the  test  section,  the  housing  was  mounted  directly  on  a  six-inch 

oil  diffusion  pump  with  its  axis,  and  the  axis  of  the  compressor  rotor, 

in  a  vertical  direction.  The  shaft  connection  between  the  rotor  anfl 

the  induction  motor  was  made  through  sn  annular  clearance  seal  O.750 

inches  in  diameter,  0.540  inches  in  length,  and  with  a  radial  clearance 

of  .002  inches.  Using  the  formula  for  a  thin  slitlike  tube,  the  molecular 

-3 

conductance  of  the  clearance  seal  is  found  to  be  1.13  x  10  liters  per 
second  for  air  at  20*C.  Thus,  if  the  pressure  on  the  motor  side  of 
the  seal  is  maintained  at  a  micron  or  less,  the  flow  from  the  region 
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of  the  motor  Is  of  the  order  of  10  micron-liters  per  second,  a 
negligible  leak  at  the  conditions  of  the  tests.  To  ensure  this 
condition,  the  motor  housing  vas  connected  to  a  two- inch  oil  diffusion  ptuqp. 
As  a  further  precaution,  a  low  vapor-pressure  silicone  grease  vas  used  as 
a  lubricant  for  the  motor  ball  bearings.  Provision  vas  made  to  admit 
gas  to  the  system  both  above  and  be  lev  the  rotor.  Seals  between  the 
system  and  the  atmosphere  were  maintained  by  means  of  welded  and 
soldered  Joints,  0- rings ,  special  gasket  seeds  for  the  ionization 
gauges,  and,  in  the  case  of  the  thermocouple  gauges  and  the  gas  inlet 
tubes,  by  red  glyptal  lacquer. 

Hie  densities  upstream  and  downstream,  of  the  compressor  rotor 

were  measured  by  hot  filament  ionization  gauges.  One  thermocouple 

gauge  vas  used  to  insure  that  the  pressure  in  the  motor  housing  vas 

below  one  micron,  and  another  vas  connected  to  the  test  section  upstream 

of  the  disk.  The  flow  rate  of  gas  admitted  into  the  system  vas  measured 

on  the  high-pressure  side  of  the  leak  valve  by  the  displacement  of 

a  column  of  mercury  in  essentially  the  same  manner  as  that  described  by 

(7) 

Dushman.  It  vas  assumed  that  the  mass  flow  rate  of  gas  through  the  leak 
valve  vas  proportional  to  the  square  of  the  pressure  on  the  high  pressure 
side  of  the  valve,  as  for  compressible  flow  at  low  Reynolds  numbers, 
and  the  flow  measurement  was  used  to  determine  this  constant  of  proportion¬ 
ality.  The  speed  of  the  rotor  vas  measured  by  a  permanent  magnet-and-coil 
tachometer  and  an  electronic  counter.  Hie  temperature  of  the  outer  race 
of  the  bearing  at  the  rotor  end  of  the  shaft  vas  monitored  with  a  thermo¬ 
couple  and  potentiometer. 
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Experlmental  Technique 

The  apparatus  was  prepared  for  testing  by  thoroughly  cleaning 
all  surfaces  to  be  exposed  to  the  vacuum  with  acetone  and  by  heating 
the  entire  apparatus  with  infrared  lamps  to  a  temperature  of  the  order 
of  100*C  for  several  hours  with  the  vacuum  pumps  operating.  At  zero 
blade  speed,  with  both  leak  valves  closed,  a  pressure  of  about  2  x  10~^mm  Hg 
was  obtained  at  the  downstream  ionization  gauge.  A  comparison  of  this 
value  with  the  reading  of  the  upstream  ionization  gauge  gave  an  app¬ 
roximate  value  for  the  flow  rate  from  leaks  and  virtual  leaks  in  the 
system,  using  the  single -blade -row  calculations  for  a  stationary  rotor 
to  determine  its  Impedance  to  gas  flow.  Since  the  density  ratio  across 
a  rotating  blade  row  is  a  function  of  the  upstream  volume  flow  rate,  the 
effect  at  these  small,  leaks  in  the  system  was  made  negligible  by  admitting 
gas  through  the  downstream  leak  valve  and  raising  the  density  level  in 
the  system.  Tests  were  conducted  at  density  levels  for  which  no  density 
drop  was  measurable  across  the  stationary  rotor. 

The  dimensions  of  the  compressor  rotor  used  are  given  in  Fig.  22. 

Angular  velocities  in  the  range  from  20,000  to  30,000  r.p.m.  gave  tip 

speeds  from  620  to  955  feet  per  second.  Ordinary  air  was  used  as  the 

test  gas  for  the  compressor  tests.  In  the  cascade  experiments ,  described 

later,  xenon  was  also  used  because  its  lower  molecular  speed  allowed 

higher  dimensionless  blade  speeds  to  be  reached.  The  pressure  dwn- 

-4 

stream  of  the  rotor  was  maintained  at  about  10  mm  Hg  by  means  of  the 
lower  leak  valve.  The  flow  of  gaa  through  the  rotor  was  controlled  by  the 
upper  leak  valve  and  the  flow  rate  was  measured  by  the  displacement  of  a 
mercury  column,  as  previously  described,  and  by  the  density  drop  across 
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the  compressor  rotor  at  zero  blade  speed. 

TVo-Rcw  Compressor 

The  two -rcw  compressor  was  built  as  a  combination  of  the  rotor 
previously  described  and  of  a  following  stator  having  a  geometry  that 
was  a  mirror  image  of  that  of  the  rotor.  Figure  24  shows  details  of  the 
arrangement  and  gives  pertinent  dimensions. 

The  spacing  between  two  rows  vas  made  as  small  as  possible  in 
order  to  minimize  the  effects  of  the  end  walls,  and  also  to  reduce  the  size 
of  the  annular  passage  designated  by  L  in  Fig.  24.  This  annulus  represents 
a  leakage  path  through  which  molecules  from  the  high  density  side  may  by¬ 
pass  the  stator  with  a  transmission  probability  of  nearly  unity.  However, 
this  throughflcw  area  was  negligibly  small  in  comparison  to  the  through- 
flow  area  of  the  stator,  and  thus  its  effect  on  the  performance  of  the 
machine  was  also  negliblble. 

Experimental  Results  with  Single-Row  Compressor 

On  Figures  25,  26,  and  27  are  displayed  the  experimental  data  and 
calculated  results  for  the  single-row  rotating  machine,  with  air  as 
the  test  gas. 

The  pressure  ratio  across  the  rotor,  for  zero  net  throughflow, 
is  shown  in  Fig.  25 •  Also  shown  are  the  experimental  data  of  Kruger 
and  Shapiro ^  &  ^),  as  veil  as  the  theoretics!  pressure  ratio  for 
zero  throughflow,  as  confuted  from  Eq.  12  smd  the  transmission 
coefficients  of  Table  1.  It  is  to  be  noted  that  the  blade  speed  S  in  the 
abscissa  of  Fig.  25  is  reckoned  at  the  surithmetic  mean  radius  between  the 
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tip  radius  of  the  blades  and  hub  radius;  correspondingly,  the  value  of 
s/b  -  .46  Is  also  reckoned  at  the  same  mean  radius.  The  present 
results  and  the  earlier  data  of  Kruger  and  Shapiro  are  in  agree¬ 
ment  to  vlthin'  the  accuracy  of  measuring  the  pressures. 

While  the  theoretical  predictions  are  in  general  agreement  with 
ti»  experimental  results,  and  while  they  certainly  follow  the  same  sort 
of  curve,  there  is  a  substantial  discrepancy,  with  the  experimental  rotor 
showing  distinctly  better  performance  than  would  be  anticipated  from 
the  theoretical  calculations.  For  instance,  at  a  blade-speed  ratio  of 
0.50,  the  theoretical  expectation  of  zero-throughflow  pressure  ratio 
is  about  2.85,  whereas  the  measured  pressure  ratio  for  zero  throughflow 
is  about  3*75* 

Fig.  26  shows,  for  several  different  values  of  the  dimensionless 
speed  ratio  &  ,  the  measured  curves  of  pressure  ratio  across  the  row 
against  volume  throughflow.  The  latter  is  expressed  in  terma  of  the 
volume  throughflow  parameter  w,  namely,  the  ratio  of  the  actual  volume 
throughflow  to  the  Volume  throughflow  incident  on  the  blade  row  from 
the  upstream  end.  The  ordinates  of  Fig.  26  have  been  normalized  with 
respect  to  the  pressure  ratio  for  zero  throughflow.  In  other  words, 
they  represent  the  pressure  ratio  for  a  given  volume  throughflow 
to  the  pressure  ratio  corresponding  to  zero  throughflow  for  the  same  <$  , 
as  displayed  on  Fig.  25.  According  to  Eq.  12,  the  curves  of  density 
ratio  vs.  volume  throughflow  should  be  straight  lines,  at  least  for 
the  limiting  case  S  CoS  ^  *  O  .  it  may  be  noted  from  Fig.  26 
that  the  experimental  data  do  indeed  follow  straight  lines  within 


the  experimental  accuracy. 

Referring  again  to  Bq.  12,  the  experimental  values  of  the 
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transmission  coefficients  ^\2  a^-<^  £g|  may  be  calculated 

from  the  intercepts  of  the  straight  line  of  pressure  ratio  vs. 
volume  flow.  Thus,  using  the  experimental  data  of  Fig.  25  and 
the  "best"  straight  lines  of  Fig.  2 6,  the  transmission  coefficients 
were  calculated  as  a  function  of  the  dimensionless  speed  ratio  . 

The  results  are  shown  in  Fig.  27.  Also  shown  on  Fig.  27  is  the 
theoretical  curve  of  transmission  coefficient  vs.  dimensionless 
speed  ratio,  as  interpolated  from  Table  1. 

The  experimental  data  in  Fig.  27  are  in  general  agreement  with 
the  theoretical  predictions.  For  positive  values  of  3  ,  the  measured 
transmission  coefficients  lie  higher  than  the  theoretical;  for  negative 
values  of  S  ,  they  lie  slightly  lower. 

All  the  experimental  results  for  the  single-row  rotating  machine 
are  essentially  summarized  in  Fig.  27,  and  it  remains  now  to  inquire 
as  to  the  possible  sources  of  discrepancy  between  the  measured  values 
of  transmission  coefficient  and  the  theoretical.  To  the  degree  that  the 
actual  conditions  of  the  experiment  duplicated  the  postulates  of  the 
theoretical  model,  we  should  expect  substantially  perfect  agreement. 
However,  there  were  several  deviations,  and  these  will  now  be  dis¬ 
cussed  and  their  significance  appraised: 

1.  The  molecular  fluxes  incident  on  the  rotor  from  the  upstream 
and  downstream  sides  were  not  necessarily  exactly  Maxwellian. 
However,  the  mode  of  operation  of  the  moving  cascade  is  such 
that  slight  departures  from  Maxwellian  conditions  should 
not  have  much  effect  on  the  transmission  coefficients. 
Accordingly,  this  is  not  felt  to  be  a  significant  source 
of  deviation  between  theory  and  experiment. 
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2.  Earlier,  In  connection  with  Eq.  1,  and  In  the  calculations  of 
2^  and  2Tg| ,  the  theoretical  calculations  vere  simplified  by 
assuming  that  Scat^i  is  approximately  zero.  This  is  not  exactly 
correct.  Estimates  of  the  amount  of  error  this  might  involve,  however, 
indicate  that  this  approximation  could  scarcely  be  the  source  of 

the  difference  between  theory  and  experiment  observed  in  Fig.  27. 

3.  In  making  the  theoretical  calculations,  it  was  assumed  that 

molecules  incident  on  the  blades  vould  be  reflected  in  a  diffuse 
manner.  That  is,  the  molecules  emitted  from  the  surface  vere 
calculated  as  though  they  came  from  a  Maxwellian  region  behind  the 
surface.  Examination  of  the  geometry  of  the  blade  rows  shows  that 
even  a  very  small  proportion  of  specular  rather  than  diffuse  reflec¬ 
tion  would  make  an  enormous  difference  in  the  transmission  coefficients. 
In  particular,  it  would  greatly  increase  the  transmission  coefficient 
^  ,  and  it  vould  greatly  decrease  the  transmission  coefficient 

’  Me  iiave  ootklog  fi™  on  which  to  judge  whether  this 
was  a  significant  factor,  but  a  comparison  of  these  results  with 
those  of  the  two-row  rotating  machine  and  with  those  for  the 
stationary  cascade  suggests  that  this  vas  not  the  essential  source 
of  discrepancy. 

4.  The  most  likely  source  of  discrepancy  between  theory  and  experiment 
in  Fig.  27,  especially  when  ve  take  into  account  the  experimental 
and  theoretical  comparisons  for  the  two-row  rotating  machine  and 

the  stationary  cascades,  is  the  fact  that  the  cascade  is  not  strictly 
two-dimensional.  Entirely  apart  from  the  variation  of  radius  and 
of  blade  speed  along  the  span  of  the  blades,  the  most  important 
factor  seems  to  be  that  in  the  actual  machine  there  are  end  wadis 
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to  the  cascade.  One  of  these  end  vails,  at  the  huh,  is 
stationary  vlth  respect  to  the  blades.  The  other  end  vail,  at 
the  casing,  Is  moving  vlth  respect  to  the  blades.  The  stationary 
end  vail  tends  to  reduce  the  values  of  both  transmission  coefficients, 
iE|g  and  ,  through  vhat  may  be  thought  of  as  a  frictional 
resistance.  If  one  examines  the  relationship  betveen  the  moving 
end  vail  and  the  blades.  It  vill  be  seen  that  any  molecules 
Incident  upon  this  moving  end  vail  vill  be  reflected  back  into 
the  blade  rev  In  such  a  vay  as  to  Improve  their  chances  of  being 
emitted  toward  the  downstream  side.  This  tends  to  Increase  the 
value  of  2^2  and  to  decrease  the  value  of  a  trend  In 

agreement  vlth  the  comparison  of  Fig.  27-  Moreover,  as  further 
Indication  that  this  may  veil  be  the  principal  explanation,  note 
that  at  a  value  of  the  dimensionless  blade  speed  ratio  S  *  0  , 
where  the  effect  of  the  moving  end  vail  vould  be  minimized,  the 
agreement  betveen  theory  and  experiment  Is  within  the  accuracy  of 
the  experimental  data.  As  further  evidence  on  this  point,  the 
effect  of  the  stationary  end  walls  at  S  =  0,  according  to  the 
Investigations  described  later  vlth  stationary  cascades  having 
various  aspect  ratios,  is  quite  small. 

experimental  Results  vlth  Two-Row  Compressor 

Figs.  26,  29,  and  30  show  the  experimental  data  and  the  eenputed 
experimental  results  for  the  tvo-rov  rotating  machine.  These  three 
figures  are  comparable  In  their  construction  to  Figs.  25,  26,  and  27 
respectively,  and  may  be  Interpreted  In  the  same  vay.  The  one  point 
requiring  clarification  has  to  do  vlth  the  theoretical  calculation  of 
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the  transmission  coefficient*  shown  In  Pig.  30.  These  were  calculated 
fra*  the  single-row  transmission  coefficients  of  Table  1  using  the  calcula¬ 
tion  procedure  described  in  ref.  5*  Vote  further  that  the  solid  curve  of 
Fig.  28  is  the  overall  zero-flow  density  ratio,  while  the  dashed  curve 
is  the  zero-flow  density  ratio  for  the  stator  alone,  found  by  combining 
the  results  of  Pig.  25  with  the  solid  curve  of  Pig.  28. 

In  assessing  the  comparison  between  the  experimental  and  theoretical 
results  for  the  two- row  rotating  sachine,  ve  may  restate  the  first  three 
consents  made  above  with  respect  to  the  similar  comparison  for  the  single- 
row  rotating  machine.  The  last  consent  made  above,  however,  requires 
modification,  inasmuch  as  the  stator  of  the  two-row  machine  had  two 
stationary  end  walls.  This  would  have  the  result  of  reducing  the  trans¬ 
mission  coefficients  and  Xg,  for  the  stator  row.  This  interpretation, 
when  applied  to  the  one-row  machine  (Figs.  25,  26,  and  27)  and  to  the 
two-row  machine  (Pigs.  28,  29  and  30),  seem*  to  make  for  a  consistent 
picture,  especially  when  it  is  realized  that  there  is  a  leakage  backflow 


around  the  stator. 
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PKEE-MOLECULE  CASCADE  TESTS 


The  third  stage  of  the  program  represented  the  most  novel 
contribution  of  the  present  research.  Its  purpose  vas  two -fold.  Ifte 
first  vas  to  enlarge  the  scope  of  existing  experimental  results  app¬ 
licable  to  the  free -molecule  compressor  by  testing  a  variety  of  different 
blade  geometries .  The  second  vas  to  test  a  new  concept  in  lew -density 
"wind  tunnels".  To  obtain  the  required  data  through  the  use  of  rotating 
blade  rows  vould  have  necessitated  the  building  and  balancing  of  a  new 
rotor  for  each  geometry  tested.  The  alternative  scheme  vas  based  on  a 
free -molecule  cascade  in  vhich  a  mean  mass  velocity  vas  inparted  to  the 
gas  on  one  side  of  a  stationary  blade  rew  by  means  of  an  uhbladed  rotor. 

Apparatus 

An  assembly  drawing  of  the  apparatus , modified  for  use  as  a  cascade 
wind-tunnel,  is  shown  in  Fig.  31a. 

The  blades  were  mounted  in  one  quadrant  of  the  cylindrical  shell 
separating  zones  I  and  ^  ,  as  shewn  In  Fig.  31b.  since  the  blades 
in  this  arrangement  are  stationary,  there  is  no  problem  of  strength 
or  balancing.  They  vere  simply  cut  out  of  thin  sheet-stock  brass 
vlth  tin  snips,  and  vere  held  in  place  in  the  blade  racks  vlth  small 
dabs  of  epoxy  resin  at  each  corner.  A  mean  mass  velocity  vas  imparted 
to  the  incident  molecules  on  one  side  of  the  cascade  by  the  high-speed 
rotor.  The  design  of  the  outer  edge  of  the  rotor,  with  two  small  ridges 
on  the  top  and  bottom,  Insured  that  all  molecules  incident  on  the  blade 
row  from  region  1  vere  emitted  from  the  outer  edge  of  the  rotor.  Therefore, 
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under  the  assumption  of  diffuse  reflection,  the  gas  in  region  3  had  a 
mean  mass  velocity  with  respect  to  the  blade  rev  equal  to  the  peripheral 
speed  of  the  rotor.  Thus  the  stream  of  molecules  emitted  from  the  rotor 
played  the  same  role  as  the  high  speed  gas  stream  in  a  rarefied-gas  vind 
tunnel. 

Procedure 

In  order  to  determine  the  transmission  coefficients,  three  sets 
of  measurements  vere  required.  The  first  was  a  calibration  of  the 
blade  row  as  a  flcw-meter  with  the  rotor  locked,  and  with  the  positive 
displacement  flowmeter  used  for  determining  Z0  ,  the  transmission 
coefficient  of  the  blade  rev  at  a  blade-speed  ratio  of  zero.  The  second 
was  a  series  of  measurements  of  the  pressure  ratio  across  the  blade  row 
at  finite  blade  speeds  and  zero  net  throughflov,  for  the  determination 
of  vs.  5.  And  the  third  was  another  set  of  runs  similar 

to  the  second,  with  the  blades  facing  in  the  opposite  direction,  to 
determine  ^  21  /  vs.  .  In  the  case  of  flat-plate  blades  with 
their  accompanying  symmetry,  it  vas  possible  to  make  the  third  set  of 
measurements  simply  by  reversing  the  direction  of  rotation  of  the 
rotor. 

In  the  original  design  of  the  apparatus,  it  had  been  intended  that 
the  resistance  to  flow  between  regions  3  and  4  and  also  between  4  and 
3  would  be  negligible  in  c oeparison  to  the  resistance  of  the  blade 
row  itself.  This  was  indeed  the  situation  in  the  latter  case,  as 
indicated  by  the  fact  that  the  pressures  in  regions  4  and  5  were  equal. 
Accordingly,  one  could  speak  meaningfully  of  a  single  downstream 
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pressure  during  the  determination  of  2fa  .  However,  the  resistance 
between  3  and  1+  was  larger  than  had  been  anticipated.  This  difficulty 
was  overcome  by  measuring  two  resistances;  first,  that  without  the 
blades  present;  and,  second,  that  with  the  blades  in  place.  From  these 
two  pieces  of  information  it  was  possible  to  evaluate  the  resistance 
of  the  blade  row  alone.  The  method  is  described  in  Appendix  B. 

The  experimental  technique  vas  essentially  the  same  as  that  des¬ 
cribed  for  the  single-stage  compressor,  with  one  addition.  As  mentioned 
earlier,  it  was  necessary  that  the  pressures  in  regions  4  and  5  be  identical 
in  order  that  the  downstream  pressure  be  a  unique  quantity.  To  this  end 
a  large  passage  was  provided  around  the  rotor  between  the  regions,  but 
an  additional  ionization  gauge  was  installed  for  region  4  and  monitored 
during  the  tests  to  Insure  that  such  equality  existed.  This  problem  of 
flow  resistance  between  the  various  regions  was,  of  course,  a  matter 
of  concern  only  in  the  calibration  test  for  determining  Z0  .  The  tests 
at  finite  blade  speed  to  determine  the  plot  of  zero-flow  density  ratio 
vs.  blade  speed  ratio  were  conducted  with  no  throughflow.  At  this  con¬ 
dition  the  resistance  to  net  throughflow  between  the  three  regions  is 
of  no  concern. 

Experimental  Results  with  Free-Molecule  Cascade 

The  experiments  performed  had  three  main  goals: 

1.  To  check  out  the  theory  for  a  variety  of  flat-plate  blade 
geometries  without  the  necessity  of  machining  new  rotors  for 
each  case. 

2.  To  investigate  the  influence  of  end -wall  effects,  in  terms  of 
the  aspect  ratio  as  an  experimental  parameter. 
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3.  To  test  non-flat-plate  blade  shapes  for  which  an  analysis 
could  not  readily  be  performed. 

All  the  experimental  data,  and  the  associated  results,  are  shown 
in  Figures  32  through  Figure  4l. 

Fig.  32  shows  the  curve  of  pressure  ratio  vs.  throughflow  for  the 
apparatus  with  no  cascade  installed.  These  data  were  employed  as 
indicated  in  Appendix  B. 

With  the  exception  of  Figure  36,  all  the  other  Figures  are  arranged 
in  three  parts,  and  each  Figure  relates  to  an  individual  blading  geometry, 
that  is,  to  a  particular  angle  oL  ,  solidity  s/b,  and  aspect  ratio  AR. 
Referring  to  Fig.  33  as  an  example,  ( =20°,  s/b  =  l/2,  AR  =  3/2), 

Fig.  33a  shows  the  curve  of  pressure  ratio  vs.  throughflow  with  the 
rotor  locked,  i.e.,  for  S  =  0.  In  conjunction  with  the  data  of  Fig.  32, 
this  led  to  the  value  of  according  to  the  method  of  Appendix  B.  Fig.  33b 

shows,  for  zero  throughflow,  the  curve  of  density  ratio  vs.  blade-speed 
ratio;  it  was  found  by  multiplying  together  the  measured  values  of 
£.2  A,  •  The  latter,  which  were  of  course  measured 

individually,  may  be  combined  vith  the  values  of  Y0  obtained  previously 
to  give  the  final  values  of  and  displayed  in  Fig.  33c. 

Effect  of  Aspect  Ratio 

Figures  33,  34  and  35  all  refer  to  the  blading  parameters  oC-  20*, 
s/b  =  l/2,  which  is  very  close  to  the  corresponding  parameters  for  the 
tests  with  the  one-rcw  and  two-row  rotating  machines.  They  differ  in 
that  they  pertain  respectively  to  aspect  ratios  of  1.5,  2.8,  and  4.0. 

Note,  by  way  of  congparlson,  that  for  the  rotating-machine  tests  the  aspect 
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rat  io  vas  about  2.4. 

The  most  interesting  feature  of  this  series  of  three  figures  is  the 
question  of  hew  the  stationary  end  vails  of  the  cascade  affect  the  trans¬ 
mission  coefficient.  For  convenience,  the  mean  experimental  curves  are 
collected  in  Fig.  36,  where  the  theoretical  curve  is  also  shown.  It  is 
evident  that  the  transmission  coefficients  for  an  aspect  ratio  of  4  lie 
nearly  at  the  asymptotic  value  for  an  aspect  ratio  of  infinity,  and  that 
even  for  an  aspect  ratio  of  2.5  the  effect  of  the  end  vails  is  not  very 
large.  The  theoretical  results  for  Infinite  aspect  ratio  are,  within 
experimental  error,  coincident  with  the  experimental  results  for  aspect 
ratio  4.  All  this  tends  to  confirm  our  earlier  appraisal  that  the  main 
reason  for  the  discrepancy  between  theory  and  experiment  in  the  rotatlng- 
machlne  tests  vas  the  pumping  effect  of  the  moving  end  wall. 

The  effect  of  aspect  ratio  observed  here  is  much  the  same  as  that 
reported  in  reference  8,  where  the  results  for  chevron  baffles  with  an 
aspect  ratio  of  5  were  nearly  coincident  with  those  for  an  aspect  ratio 
of  infinity. 

Effect  of  Spacing-Chord  Ratio 

Figures  35,  37,  38,  and  39  are  all  for  a  blading  angle  cJ.  =  20* 
and  with  an  aspect  ratio  large  enough  (  >4)  to  be  considered  nearly 
infinite.  They  differ  in  that  each  is  for  a  different  spacing-chord 
ratio,  s/b ,  covering  the  range  from  1/4  to  3/2 .  In  every  case  the 
measured  transmission  coefficients  are  in  close  agreement  vlth  the 


theoretical  values 
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Effect  of  Blade  Curvature  ' 

Figures  40  and  Ul  shew  experimental  data  for  cambered  blading,  for 
which  no  theoretical  results  are  available.  Each  of  the  two  cambered 
design  is  based  on  a  straight-chord  line  with  oi  =  20°,  s/b  =  l/2, 
and  the  results  are  therefore  to  be  compared  with  those  in  Fig.  35* 

With  the  camber  of  Fig.  UO,  both  the  zero-flcw  pressure  ratio 
and  the  volume -flew  capacity  are  less  than  with  straight  blades  having 
the  same  OL  and  s/b.  Therefore  this  direction  of  camber  is  undesirable, 
at  least  in  the  range  of  oC ,  s/b,  and  camber  of  the  tests. 

With  the  opposite  camber,  Fig.  Ul,  the  zero-flcw  pressure  ratio 
is  slightly  larger  than  that  for  flat  plates,  but  the  through-flcw 
capacity  is  only  about  half  as  large.  Therefore  this  direction  and  amount 
of  camber  also  seems  to  offer  no  substantial  advantage. 
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APPENDDC  A 

DERIVATION  OF  FUNCTIONS  USED  IN  ANALYTICAL  SOLUTION 

Initial  Incidence  Distributions 

As  explained  earlier,  the  analytical  solution  may  be  broken  up 
into  two  major  steps,  the  first  of  which  is  the  calculation  of  the 
so-called  initial  incidence  distribution.  By  reference  to  Fig.  A-l,  it 
is  seen  that  molecules  Inpinging  upon  the  blade  row  from  region  1  may 
suffer  one  of  three  fates.  They  will  either  pass  directly  through  the 
blade  row  without  a  collision  with  the  surfaces  and  thus  be  of  no  further 
interest,  or  else  they  will  collide  with  one  of  the  two  surfaces.  The 
items  of  importance  in  this  Btage  of  the  calculation  are:  first,  the 
fraction  which  pass  through  the  row  directly,  which  we  call  ; 
and  second,  the  fractions  which  collide  with  areas  dx  and  <ky  on  the  blade 
surfaces  at  varying  X  and  y  along  the  blades .  These  distribution 


Fig.  A-l 


This  problem  will  be  approached,  la  the  notation  of  Fig.  A-2 , 


Fig.  A-2 

by  first  calculating  the  number  of  molecules  incident  upon  area  <Lx  at  x 
which  enter  the  blade  row  by  crossing  the  area  <£?  at  ?  ,  and  then  Integrating 
across  the  opening  between  the  blades  from  2  =  0  to  2  =  s  .  Since  the 
Incoming  molecules  are  uniformly  distributed  in  space,  every  d.2  has  the 
same  molecular  flux.  Therefore  the  problem  is  to  determine  for  a  given 
blade  speed,  what  fraction  of  molecules  crossing  <Ai  are  contained  in  the 
plane  angle  d.Q  inclined  at  the  angle  B  to  the  normal.  A  knowledge  of 
the  blade  geometry  will  then  relate  the  x,  2,  and,  B  ,  thus  enabling 
one,  at  least  in  principle,  to  perform  the  necessary  integrations. 

From  the  point  of  view  of  an  observer  on  the  blade  row,  the  molecules 
Incident  upon  the  row  from  region  1  will  have  a  Maxwellian  velocity  dis¬ 
tribution  superposed  on  a  mean  mass  velocity  vector,  specified  by  V  and 
jl  in  Fig.  A— 3 •  Let  U  represent  the  molecular  velocity  component 
along  |  ,  IT  along  ij  ,  and  W  along  (normal  to  paper). 
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© 


Fig-  A-3 


Hence  the  number  of  molecules  per  unit  volume  having  velocity  components 


lying  between  XL  and  XL  +■  ,  V  and  V'Mlir^and  vr  and  ter i-dw  is 

.  v  (9) 

given  by 


+  (v-Va*fl)  4-  J\ 


oLu  drfdw 


(Eq.  A-l) 


Transforming  from  cartesian  velocity  coordinates  U,7VjUT  to  cylindrical 
coordinates  v,b,  and  ur  ^cording  to  the  relations 

U  =  U  ars  0  (Eq.  A-2a) 


VS  U  sUl  0 


(Eq.  A- 2b) 


1A1=  vr 


(Eq.  A-2c) 


the  volume  in  velocity  apace  dll  dUrdu;  becomes  UcBDdfidlUr  and 
Eq.  A-l  becomes 


Uduaedur  «*»•  A-3> 
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Now,  ail  the  molecules  which  cross  a  surface  per  unit  area  per  unit  time 
with  velocities  between  U  and  U+dU  in  the  angle  d&  at  angle  9  to  the 
normal  must  have  come  from  an  Imaginary  cylinder  behind  the  surface  with 
slant  height  and  altitude  U  COS  Q  .  Thus  the  number  of 

these  molecules  crossing  unit  surface  per  unit  time  is  given  by 


(lTT?T) 


a 


-  6  -Vcosy&f-v  e  -Vs^  /if+w2]^ ; 


‘^eaUdodw  (Eq-  A-4> 


If  this  is  now  integrated  with  respect  to  -ui  from  -  oo  to  -t-oo  .and  with 
respect  to  U  from  0  to  oo  }ve  obtain  the  total  number  of  molecules  crossing 
the  surface  per  unit  time  and  per  unit  area  with  velocity  vectors  lying 
in  dd  at  d  : 

'n.tfzirr 


2>r 

where 


$  =  V//ZST 
$,  =  $  oos  (e-/i) 
$2  =  3^  (0-/i) 


C05  0  cLO  (Eq-  a- 5) 


(Eq.  A-6) 


Dividing  now  by  the  total  flux  of  molecules  crossing  the  plane  to  get  the 
fraction  of  the  total  with  velocities  lying  in  d&  at  9  ,  one  obtains 


cos  9  (Eq-  A-7) 


i  “s 

»e  + 

(4+s.£)e 

c 

fF  3orsA| 

& 
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If  now,  both  the  blade  surfaces  and  the  opening  between  the  blades 
are  broken  up  Into  a  finite  number  of  sections  as  shown  In  Fig.  A-U,  the 
angle  subtended  by  the  £th  small  area  ix  from  a  point  z  In  the  center 
of  the  ith  interval  is  denoted  by 


Ad 


Where,  in  general 


hi 


=  6;  ;  -  6, 


Fig.  A-U 


0  =  4"&n. 


-I 


(s/b)(i/s)  +  (X  A)  Cosot 
(x/b)  ^  °L 


(Eq.  A-8) 


Using  Eqs.  A-8,  Eq.  A-7  may  be  integrated  between  0^  j  and 
for  a  fixed  position  l  on  the  scale  of  Z  .  One  may  then  sum  over  all 
the  t  sections  of  the  opening  between  the  blades  in  order  to  obtain  the 
total  fraction  of  molecules  from  the  region  considered  which  are  directly 
incident  upon  dx^ .  Having  done  this  for  all  dxi  and  dy  i ,  one  has  a 
stepwise  approximation  to  i/  and  l/.  .  The  fraction  transmitted 

X  7 

directly  from  region  1  to  region  2  is  calculated  in  a  similar  way. 
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Ultimate  Transmission  Probability  of  Reflected  Molecules 


In  performing  the  second  stage  of  analysis,  that  of  determining  the 
ultimate  probability  of  a  molecule's  getting  through  to  region  2  after 
being  emitted  for  a  known  position  along  the  blades,  it  will  be  remembered 
that,  in  general,  it  was  necessary  to  solve  two  simultaneous  Integral 
equations  for  along  with  <^z  which  were  of  the  form. 


crxz  -  Cxl  iffy 

and 

F(% *■)<** 


(Eq.  A-9&) 


(Eq.  A- 9b) 


where  was  defined  as  the  fraction  of  molecules  emitted  from 

dx.  at  X  which  impinge  directly  upon  the  area  dy  at  y  ,  and  Cxtthe 
fraction  of  molecules  emitted  from  dx  at  X  which  reach  region  2  directly 
without  any  further  collision  with  the  blade  surfaces.  and 


are  defined  similarly. 

For  the  case  under  consideration,  certain  simplifications  are  possible. 
First,  if  the  cosine  law  of  reflection  is  valid >  F(a/!l)  =  F  (%  x)  •and 
second,  for  the  case  of  flat  plate  blading,  symmetry  enables  one  to  say, 
for  equal  values  of  x  and  y  ,  rxz  ~  ,  and,  similarly^ 

Furthermore,  since 


I 


+ 


(Eq.  A- 10) 


it  is  possible  to  reduce  the  two  simultaneous  equations  to  a  single  one 
in  the  form, 


<r  * 

xt 


l- 


^1  F  (*<!))  ^ 


(Eq.  A-ll) 
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Slnce  C£  (x)  and  (y)  are  the  same  function,  this  equation  may  be 
solved  by  iterative  methods  once  the  functions  for  F(x,  y)  and  Cx, 
are  known. 

For  the  calculation  of  Cx|,  refer  to  Fig.  A-5. 


Fig.  A-5 

Assuming  the  cosine  law  of  diffuse  reflection,  the  fraction  of  molecules 
emitted  in  an  angle  d&  Inclined  at  an  angle  9  to  the  normal  is 

j  Cos  9  c(9 

Therefore  the  fraction  emitted  in  the  angle  "/  is 


r 


72 


C  =  i 

Wi  l 


00$ 


s  diQ  =  i  ( i-coi  y) 


(Eq.  A- 12) 


1./ 

1  ' 


From  the  geometry  of  Fig.  A-5, 


/  = 


-I 


+  T 


<LX)% 


‘/2 


(Eq.  A-13) 
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Accordingly, 


Tor  the  calculation  of 


refer  to  Pig.  A-6. 


Again,  by  the  assumption  of  diffuse  reflection, 

w  -j  0  dO 

Referring  to  Fig.  A-6  for  the  geometrical  configuration, 


Also, 


filft  = 


cel  8 

r 


(Eq.  A-14) 


Is  given  by 


(Eq.  A-15) 
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Hence 


i  (if  4i) 

F£x,^«  - — —  (Kq.  A-l£) 

b 
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APPEHDIX  B 

DETTOMIHATIOW  OF  IN  CASCADE  TESTS 

As  described  earlier,  the  free-molecule  cascade  with  zero  throughflow 
enables  one  to  measure  directly  z/^-o  “d  /?0  as  a  function  of 
blade-speed  ratio  £  .  Therefore,  In  order  to  obtain  the  magnitude  of 
the  transmission  coefficients  T|$  and  ,  one  must  measure  ^  for  the 
stationary  cascade.  This  is  done  by  measuring  pressure  ratio  vs.  volume 
throughflov,  and  using  the  relationship  "n^/n.,  *  I  -  (yj/'Z0')  . 

In  the  cue  of  the  rotatlng-rov  machine,  this  was  a  straightforward 
task  since  the  only  appreciable  resistance  to  flow  between  the  upstream 
and  downstream  pressure  gauges  was  the  blade  row  itself.  In  the  case  of  the 
stationary  cascade,  such  was  not  the  case.  Hcwever,  the  determination  of 
was  arrived  at  through  the  following  considerations.  Fig.  B-l  is  a 
diagram  of  the  stationary  cascade,  showing  the  several  resistances  to  flow 
between  the  two  pressure  gauges. 


Fig.  B-l 
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In  general,  if  one  has  several  flow  resistances  in  series  between  two 
points,  it  is  possible,  from  a  measurement  of  the  flow  at  a  given  pressure 
ratio  between  the  two  points,  to  determine  the  value  of  any  individual 
resistance  if  all  the  others  are  known.  What  is  desired,  in  this  instance, 
is  the  value  of  the  resistance  from  2  to  3  when  the  blades  are  in  place. 
This  requires  a  knowledge  of  the  resistances  from  1  to  2,  3  to  4,  and 
4  to  5*  Of  these,  only  the  resistance  from  3  to  4  was  not  known  and  not 
readily  calculable.  Therefore  it  was  determined  by  means  of  a  preliminary 
experiment. 

The  general  equation  for  the  flow  from  a  region  i  to  a  region  J  may 
be  written  as 


Wij  -  Wf  - 


(B-l) 


where  M  signifies  molecular  flux;  the  transmission  coefficient  from 
1  to  J  ;  and  the  net  throughflux  ratio  based  on  the  Incident  flux 
HA;  for  region  i . 

In  the  absence  of  moving  parts,  it  is  further  required  that 


Al  1ii  ' 


(B-2) 


Applying  these  equations  to  the  successive  stations  of  Fig.  B-l,  we  get 


(B-3a) 

Us/'n1  -  1- 

Ia4 

(B-3b) 

/r.3  *  1 

-  w3  /r» 

(B-3c) 

-  Wn  /t,y 

(B-3d) 
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Slnce  the  mass  flow  is  constant  along  the  flow  path, 


W;*;  “  K 


(B-4) 


where  K  is  a  constant.  Now,  since  the  positive  displacement  flew  meter, 
together  with  the  pressure  gauge  at  1,  yields  the  value  W. ,  It  is  convenient 
to  eliminate  K  by  writing 


U,  Ai 

t  X 


U  l 


(B-5) 


Since  /Aj  ve  may  approximately  set  ^  =  0  .  Also, 

since  £  1  jit  follows  that  iS  1 

By  approximating  the  flew  between  stations  1  &  2  and  between  2  &  3 

as  flow  between  parallel  plates,  the  values  It,*  =  0.13  and  ^25“  60 

(7) 

are  obtainable  from  the  literature  .  Since  these  passages  represent 
relatively  small  resistances  to  the  flow  as  compared  to  the  passage  from 
3  to  4,  small  inaccuracies  in  these  transmission  coefficients  have  a  very 
slight  effect  on  the  determination  of  the  resistance  from  3  to  4. 

With  these  quantities  known,  a  measurement  of  at  any  given  value 
of  Hj/aj  is  sufficient  to  solve  for  using  the  following  combination  of 
Equations  B-3  and  B-5: 

* 

*5- /  wat,  %  a  ,  wj_ a,/ ,  hYi.  \ 

Numerically^  A,  /Atw  A,/As *  W7 .  from  Fig.  32,  at  K5  /a(  *  0.  5?  W,  =  0. 0  2<| . 
Thus  one  may  finally  calculate  ~  0.  35 \ . 

With  this  value  of  tests  may  be  run  with  the  blades  in  place,  and 
application  of  the  equation  above  to  the  data  furnishes  values  of  Z0 
(i.e. ,  Z^)  for  each  blade  geometry. 
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TABLE  I 


THEORETICAL  VALUES  OF  TRANSMISSION  COEFFICIENT,  X 


Note:  Positive  values  of  S  correspond  to  X  ni  negative  values 
correspond  to  Z  21* 


(a)  s/b-1/4 


(b)  s /b*l/2 
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TABLE  I  (continued) 


THEORETICAL  VALUES  OF  TRANSMISSION  COEFFICIENT, 2 


(c)  s/b-1 


a 


s 

!  io8 

i  20° 

30°  i 

!  45°  ] 

1  6 

•  5097 


0 

•  4773 


1.00 

-1.00 

•  3311 

•  0446 

•  5542 

•  0968 

2.00 

-2.00 

•  5258 

•  0234 

3.00 

-3.00 

■Bjm 

nm 

BnTH 

.8580 

•0405 

amisi 

•  8433 
.0135 

•  8365 
.0376 

.7048 
•  0703 


.3569 


.7143 
•  2962 
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.2497 
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.1462 
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•  4031 
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•  3498 
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TilHORLiT  I  CAL  VALUES  OF  TRANSMISSION  COEFFICIENT,  S 


(e)  s/b*2 


a 

s 

10° 

20° 

O 
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TO 
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0.00 
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FIG.  2 


FIG.  3 

THEORETICAL  TRANSMISSION  COEFFICIENTS  FOR  SINGLE 
ROW,  s/b  =  1/4 


FIG.  4 

THEORETICAL  TRANSMISSION  COEFFICIENTS  FOR  SINGLE 
ROW,  s/b  =  1/2 
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FIG.  6 

THEORETICAL  TRANSMISSION  COEFFICIENTS  FOR  SINGLE 
ROW,  s/b  =  3/2 


S  =  VB  A/2RT 
FIG.  7 


THEORETICAL  TRANSMISSION  COEFFICIENTS  FOR  SINGLE 
ROW,  s/b  =  2 


S  =  VbV2RT 
FIG.  9 


THEORETICAL  VALUES  OF  W  FOR  ZERO  PRESSURE  RISE 
(Z,2-Z2I),  AND  OF  n2/n,  FOR  ZERO  FLOW  (Z,2 /Z2I),  FOR 
s/b  a  1/2 


THEORETICAL  VALUES  OF  W  FOR  ZERO  PRESSURE  RISE 
(ZI2-Z2|),  AND  OF  n2/n,  FOR  ZERO  FLOW  (ZI2/Z2I),  FOR 

s/b  =  I 


THEORETICAL  VALUES  OF  W  FOR  ZERO  PRESSURE  RISE 
(  2i2-S2|  ),  AND  OF  n2/n,  FOR  ZERO  FLOW  (ZI2/S2I),  FOR 


FIG.  12 

THEORETICAL  VALUES  OF  W  FOR  ZERO  PRESSURE  RISE 
(£i2-Z2i),  AND  OF  n2/n,  FOR  ZERO  FLOW  (Sl2  /Z2,),FOR 
s/b  =  2 
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FIGURE  19 


FIG.  20 
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STATIONARY  CASCADE  RING 
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